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Parity measurement is a central tool to many quantum information processing tasks. In this 
Letter, we propose a method to directly measure two- and four-qubit parity with low overhead 
in hard- and software, while remaining robust to experimental imperfections. Our scheme relies 
on dispersive qubit-cavity coupling and photon counting that is sensitive only to intensity; both 
ingredients are widely realized in many different quantum computing modalities. For a leading 
technology in quantum computing, superconducting integrated circuits, we analyze the measurement 
contrast and the back action of the scheme and show that this measurement comes close enough to 
an ideal parity measurement to be applicable to quantum error correction. 


The ability to measure a quantum system in a high 
fidelity and quantum non-demolition (QND) way is fun¬ 
damental to most aspects of quantum information pro¬ 
cessing (QIP). In the circuit quantum electrodynam¬ 
ics (cQED) community, great success has been achieved 
in qubit readout by linear amplification and homodyne 
detection of the signal in a dispersively coupled mi¬ 
crowave resonator (cavity) [1-4]. However, the macro¬ 
scopic size and complicated circuitry required for this 
readout scheme is a major obstacle to scalability. For 
this reason a simpler, scalable, high-fidelity, QND single¬ 
qubit readout scheme based on threshold photodetection 
was recently introduced [5]. 

In addition, QND readout of multi-qubit operators is 
increasingly important in contemporary QIP. In particu¬ 
lar, readout of the parity of multiple qubits has applica¬ 
tions to quantum error correction [6] such as the surface 
code [7], to quantum phase estimation [8], to the imple¬ 
mentation of multi-qubit gates [9, 10], and to entangle¬ 
ment generation [11-13]. Parity measurement of two su¬ 
perconducting qubits has been demonstrated ]12, 14], as 
has Bell-state measurement ]15] and parity measurement 
of a cavity state using a superconducting qubit ]16]. How¬ 
ever, adapting these protocols to more than two qubits 
requires a significant increase in circuit complexity, ei¬ 
ther in the number of microwave resonators necessary for 
amplification-based direct parity measurement ]17, 18], 
or in the number of costly entangling gates in gate-based 
parity measurement. 

In this Letter, we propose a protocol for QND readout 
of the parity of multiple qubits coupled to one resonator, 
over a timescale comparable to that of a single entangling 
gate, following Ref. ]5]. We describe the physical model, 
briefly review Ref. ]5], and show how the procedure gen¬ 
eralizes to an fV-qubit single-shot parity measurement. 
Finally, we analyze the main sources of error. 

We describe parity measurement on N qubits, each 
coupled to the same single-mode cavity, which is also 
coupled to a tunable photon counter. Each qubit Q„ 
of transition frequency loq^ couples to the cavity via a 


Jaynes-Gummings interaction in the dispersive regime, 
leading to the full cavity-qubit Hamiltonian (h = 1) 

H = Hd+ujco^ a+^ f XQ„a'l'a- ^ ^ 

n ^ ^ 

where — ojq„) is the dispersive shift for 

coupling strength gq^, and Hd = A{t) (a -I- a'^) is a time- 
dependent classical drive that controls the cavity. For 
this Hamiltonian the computational basis states are de¬ 
fined as the eigenvectors of J2n assume that 

the qubits are far enough detuned that cavity-mediated 
qubit-qubit coupling is negligible. This model is suffi¬ 
ciently general to encompass many qubit architectures 
]1, 3, 12, 17, 19, 20]. 

Our protocol can be applied to any physical system 
involving a threshold photon counter and strong disper¬ 
sive coupling between qubits and electromagnetic radia¬ 
tion. In the microwave regime, where photon counting 
is possible with a Josephson Photomultiplier (JPM) ]21], 
other examples include NV centers in diamond ]22], Ry¬ 
dberg atoms ]23, 24], and lateral quantum dots ]25]. At 
higher frequencies, where photon counters are common¬ 
place, candidate systems include trapped atoms ]26] and 
self-assembled quantum dots ]27]. 

As described in Ref. ]5], the cavity-mediated QND 
measurement of a single qubit with a photon counter is 
performed in a three-stage protocol. First the cavity, 
initialized to vacuum, is driven at its shifted resonance 
frequency for the qubit in its excited state, populating 
the cavity conditionally upon the state of the qubit. In 
the following measurement stage, the photon counter dis¬ 
tinguishes between the two conditional cavity states (one 
of which is above and one below its threshold), and there¬ 
fore detects the state of the qubit. Finally, an additional 
microwave drive resets the cavity approximately to the 
vacuum state for further computation. Two properties 
of the counter are crucial: it should be a square-law de¬ 
tector that responds to the total energy only (insensitive 
to phase), and it should have threshold behavior at a 
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suitable number of photons that is well between vacuum 
and the selectively excited state of the cavity. Number 
resolution and single-photon sensitivity are not required. 


Microwave 
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single-qubit readout these are projectors onto the qubit 
ground and excited states, while for more than one qubit 
they are sums of projectors onto a set of states spanning 
each parity subspace. 

Now consider two qubits coupled to the cavity with 
equal dispersive shifts. For the computational basis 
states the total dispersive shift on the cavity is xq = 
±2xq when both qubits are in the same eigenstate (even 
parity), while XQ = 0 when the qubits are in differ¬ 
ent eigenstates (odd parity), as illustrated in FIG. 1(b). 
By driving with frequency wd = wc, we entangle the 
qubit parity subspaces with distinguishable cavity states 
(as was done for single-qubit readout), implementing 
the unitary of Eq. (17). At the end of the drive pulse 
A{t) = Oq cos (wd^ + </')0(iD — t), the qubit-dependent 
cavity occupations are 

I |2 f ao V 1 - cos (ADto) I |2 /ao 

=(^) -^-’ i“°i “It 

where Aq = Cjc ~ <^d- Setting wq = wc gives |Ad| = 
2xq for even-parity states, and therefore |aEp = 0 at 
= t/xq (^D = 0 for the odd states). As = 0, the 
even-parity states are indistinguishable, and thus, while 
inter-subspace coherence is reduced to (ao| Qie), intra¬ 
subspace coherence is protected, unlike in direct parity 
measurement [12]. 

For N > 2 qubits, the degeneracy within the odd and 
even subspaces splits and we can no longer perform par¬ 
ity measurement by a single frequency cavity drive. How¬ 
ever, with 



FIG. 1: (a) Schematic of the experimental design showing the 
cavity, N coupled qubits, and the photon counter (a JPM in 
this case). (b) and (c) Illustration of the dispersive shifts on 
the cavity for two- and four-qubit computational basis states. 

During the drive stage, the reduced Hamiltonian of the 
cavity coupled to a single qubit. 

He = ojcd^d + A{t) (a + a'^) , (2) 

is that of a single mode oscillator with resonance fre¬ 
quency oje = ujc + Xq, where xq = sxq encodes the 
state of the qubit through s = ±1, with s the respective 
eigenvalue of dj,. By applying a carefully timed classical 
drive at frequency wd = wq ~ XQ, we put the cavity in a 
high-amplitude coherent state if the qubit is in its excited 
state, but in the vacuum state if the qubit is in its ground 
state, thus dephasing the qubit. The cavity is populated 
conditionally upon the qubit having odd parity (excited 
state), effectively applying the unitary 

UB = P^<E>lc+P§^DiP), (3) 

which transfers qubit parity information to the cavity, 
where D{f3) is the displacement operator. Here Pq^^ are 
the projectors onto the even/odd parity subspaces. For 


A(t) = tto ^ cos (wdT + 0) 0<t<tn (4) 

i 

where wd; = wq -b XQ; ^re the dispersively-shifted cavity 
frequencies for each band of the odd-parity subspace, we 
apply the unitary of Eq. (17) by simultaneously driving 
all odd-parity spectral lines resonantly with a multi-tone 
drive (see the supplementary material for further details). 

For four qubits, basis states with odd parity pro¬ 
duce dispersive shifts of ±2 xq (blue and red odd-parity 
bands), while basis states with even parity cause a dis¬ 
persive shift of 0 or ±4xq, as shown in FIG. 1(c). There¬ 
fore, with the two drive frequencies wd, 2 = i 2 xq we 
can simultaneously drive both odd-parity spectral lines. 
The cavity occupations for the four-qubit parity bands 
are shown in FIG. 2 (see the supplement for their an¬ 
alytic form). Grucially, |aEp = 0 for all even-parity 
states at to = tt/xq, while [aop is equal for all odd- 
parity states. Therefore, qubit parity information has 
been mapped onto cavity occupation in a way that does 
not discriminate between states of the same parity, while 
reducing inter-subspace coherence to (ao| Q^e)- 

In the measurement stage we distinguish between the 
qubit-parity dependent cavity states in a frequency in¬ 
sensitive way (to avoid intra-subspace decoherence) by 
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FIG. 2: Cavity photon number as a function of the length of 
the drive pulse, to, for the four qubit parity bands. Xq/’’^ ~ 
10 MHz, so the optimal drive time is to = 100 ns. 

using a high bandwidth photon counter. We tune the 
counter on resonance with the bare cavity frequency ojc- 
For two qubits, the counter and cavity are resonant if the 
qubits have odd parity, maximizing the detection prob¬ 
ability of |ao)- For four qubits, the counter is resonant 
with the cavity if the qubits are in the center even-parity 
band; however, this is not an issue because the cavity is 
unoccupied for even-parity qubits. The counter will be 
± 2 xq off resonance from the cavity if the qubits have 
odd parity, symmetrically between the two odd-parity 
cavity frequencies. For a counter with bandwidth > dyq 
this ensures the detection probabilities for all odd-parity 
states are nearly identical, minimizing intra-subspace de¬ 
coherence. 

For = 4, the bright count detection probability 
for the two odd-parity bands is shown in FIG. 3, as is 
the false “even-state detection probability” due to dark 
counts. As expected, the two bright count curves are 
identical. These curves are calculated by solving the mas¬ 
ter equation with relevant experimental parameters (for 
further information see the supplements), where, as an 
example of a threshold photon counter we have chosen 
the JPM, which when coupled to a cavity approaches 
unit detection probability with large bandwidth (;g> yq) 
[28], and is a true dichotomic detector that distinguishes 
between the vacuum and its complement [29, 30[. Similar 
detectors exist in the optical [31, 32[ and near infrared 
[33[ regimes. 

As a figure of merit we consider the parity measure¬ 
ment contrast, defined as 

C {tu) = P (|aop, tyi) — P (|aEp, ^m) , (5) 

the difference in detection probability between odd/even 
parity states, which is a function of the measurement 
time tu- The drive time tu is chosen such that «£ = 0, 
so the parity measurement contrast reduces to the detec¬ 
tion probability of the odd-parity cavity state minus the 
dark count probability. In reality there will be thermal 
photons in the cavity. However, because of the threshold 



FIG. 3: Bright count probability for the four-qubit odd-parity 
bands, dark count probability for even parity, and measure¬ 
ment contrast as functions of measurement time. 

nature of the photon counter it is not necessary to have 
an initially empty cavity to maximize measurement con¬ 
trast; it is sufficient that the detection probability has a 
sharp threshold between aq and ag- This is important 
for the JPM, where, while in principle the threshold can 
be set at zero occupation [28, 29[, this is often difficult to 
engineer. Nonetheless, a suitable threshold between aq 
and aE is easily obtained [5[. 

The parity measurement contrast is limited by 
misidentification of the parity, which occurs either due to 
a dark count or the nonzero probability of not detecting 
the state [aq). These errors are controlled by the bright 
count/dark count ratio of the photon counter and the 
cavity occupation jaop, and control over these param¬ 
eters is sufficient to obtain contrast arbitrarily close to 
unity. As seen in FIG. 3, measurement contrast that ap¬ 
proaches 95% is achievable, compatible with the readout 
threshold for error correction [7[[37[, in a 40 ns time- 
frame with experimentally relevant bright count/dark 
count ratio and jaqp. 

In the reset stage we again implement the qubit- 
dependent cavity displacement of Eq. (17), but with 
/3 = —aM- If the qubits have odd parity, the cavity 
begins approximately in the coherent state jaivi) and is 
returned to vacuum, where |aM| can be calculated given 
the input magnitude jaoj and the form of the detection 
back action on the cavity [30[. For a JPM, after mea¬ 
surement the cavity is not a coherent state, and reset 
is imperfect. However, single-shot parity measurement 
is unaffected by reset error, and though repeated mea¬ 
surements can be affected, for experimentally relevant 
parameters reset error is on the order of 1% [5[. 

In a realistic experiment it is unlikely that all qubit- 
cavity dispersive shifts will be identical. Therefore, we 
examine the robustness of our protocol under variations 
in the dispersive shifts, defining the dispersive shift er¬ 
ror e, such that xQi = XQ> and yq^ = yq -|- e for two 
qubits. For four qubits there are three dispersive shift 
errors: £ 2 , £ 3 , and £ 4 . We assume this error is small. 
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such that Ci/xQ <C 1. For a superconducting qubit with 
coherence time 5 /iS, the uncertainty in XQ due to homo¬ 
geneous broadening is no more than ± 0 . 1 %, and thus 
ki/xql < 0.2%. Mismatch in xq causes a reduction 
in measurement contrast and intra-subspace decoherence 
(dephasing of a superposition of qubit states of the same 
parity), and we provide quantitative estimates of these 
effects, leaving further details to the supplement. Errors 
in the drive frequencies wd, cause similar errors as xq 
mismatch, but can be corrected for during tune-up. 

Imperfect resonance between the applied drives and 
the cavity due to xq mismatch leads to a reduced |aoP 
and an increased |aEp- For either two or four qubits, 
the reduction in measurement contrast is second order 
in ei/xqj commensurate with a square-law detector. As 
shown in FIG. 4, even for 10% xq mismatch the resulting 
measurement error (even-state detection probability) is 
on the order of 1 %. 

The intra-subspace decoherence caused by xq mis¬ 
match is due to qubit states with the same parity being 
entangled to cavity states with different phases, mak¬ 
ing them distinguishable (in principle). The resulting 
decoherence is at most second order in ei/xQ- Shown 
in FIG. 4 is the coherence of odd- and even-parity two- 
qubit superposition states, as quantified by the relevant 
off-diagonal matrix element of the reduced two-qubit den¬ 
sity matrix. Though the decoherence of an odd-parity su¬ 
perposition is nontrivial for larger values of e/xq, perfect 
reset returns full coherence by depopulating the cavity in 
a phase insensitive way. 

However, any cavity decay prior to perfect reset will 
cause irreversible coherence loss, as would the decay of 
residual photons after imperfect reset. Up to the limit 
of strong xq-mismatch (e approaching xq)i intra¬ 
subspace coherence decays to exp(—A^c(l — cos(7re/xq))) 
in the steady state, where Nq is the average photon 
number in the cavity prior to decay, hence requiring 
e < Xq/V^C- For a sufficiently high-Q cavity, decay 
during the short measurement time is unlikely, and post¬ 
reset cavity decay is a more significant source of loss. For 
the parameters considered here, the residual photons af¬ 
ter imperfect reset result in coherence loss that is less 
than 1% in the worst case, and is therefore inconsequen¬ 
tial. 

For two qubits with xq mismatch, no intra-subspace 
decoherence occurs during the measurement stage as cav¬ 
ity states corresponding to qubit states of the same parity 
have equal detection probability. For four qubits with xq 
mismatch the detection probabilities for states within the 
same parity subspace can be different, with their differ¬ 
ence scaling as (ef — e|)/xQ. This changes the basis of 
measurement, which changes the output qubit state at 
the end of the protocol. This effect can be mitigated by 
a square-law photon counter that quickly saturates with 
photon number above threshold, such as the JPM, for 
which the error in the output qubit state is negligible. 



FIG. 4: Even-parity detection probability, odd-parity coher¬ 
ence, and even-parity coherence as functions of the XQ mis¬ 
match for two-qubit parity measurement. e/xQ ranges from 
1% to 20%, well within experimental expectations. xq/Stt = 
5 MHz, and the drive power oq is such that |ao|^ = 9. 


Higher order effects beyond the dispersive Hamiltonian 
can affect the parity readout protocol presented here. As 
a result of the formation of cavity-qubit dressed states, 
there will be residual cavity occupation for even qubit 
states, which reduces parity measurement contrast. For 
four qubits the contrast is reduced to « 92% (see supple¬ 
ment). However, contrast can be improved by increasing 
cavity-qubit detuning, and/or using better drive pulse 
sequences. The protocol’s QNDness will also be affected 
by higher order terms; however, as shown for single-qubit 
readout in [5], this change in QNDness is minimal. As 
the measurement stages of single-qubit readout and par¬ 
ity readout are effectively identical, this is also true for 
parity readout. 

Our protocol involves an entangling operation that 
maps qubit information to the cavity such that destruc¬ 
tive readout of the cavity non-destructively determines 
the qubits’ parity. It provides a tailored and efficient 
quantum circuit that maps parity information onto an 
ancilla (the cavity), and then measures the ancilla in a 
way that is insensitive to qubit-resolving information. It 
occupies a middle ground between the direct, amplifier- 
based parity measurement of [17, 18] and gate-based par¬ 
ity measurement. The advantage of our protocol over 
gate-based protocols is the reduced number of entangling 
gates required: one for our proposal versus four for the 
traditional gate-based four-qubit protocol. Our protocol 
also requires fewer cavities than direct four-qubit parity 
readout [17, 18[. 

In conclusion we have presented a high-fidelity, scal¬ 
able, QND protocol for parity readout of two or four 
qubits via photon counting of a dispersively coupled cav¬ 
ity. Measurement contrast is limited by the bright/dark 
count ratio of the photon counter, and with the limited 
optimization in this Letter approaches 95%. Our proto¬ 
col introduces no decoherence of qubit states with the 
same parity, and is robust against the major sources or 
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error. Future work will focus on protocol optimization 
and higher order corrections. 
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SUPPLEMENTARY MATERIAL 

In the following supplemental sections we expand on 
the details of the calculations and simulations presented 
in the main text. While not essential to understand¬ 
ing the results presented in the main text, these details 
complete the explanation of the results in a more math¬ 
ematically rigorous way. 

In section we present the details of the analytic deriva¬ 
tion of FIG. 2 of the main text, which shows the driven 
cavity occupation as a function of time, and how this is 
used to create a qubit-parity dependent cavity popula¬ 
tion. In section we describe the master equation formal¬ 
ism used to simulate the interaction of the cavity with 
our example threshold square-law photon counter, the 
Josephson Photomultiplier (JPM). In section we expand 
on the error analysis of the main text to describe the ef¬ 
fect of XQ mismatch. In section we discuss the effect of 
this XQ mismatch when cavity photon loss occurs. Fi¬ 
nally, in section we describe how consideration of the 
full Jaynes-Cummings Hamiltonian affects the results of 
the main text. 


DERIVATION OF THE QUBIT-PARITY 
DEPENDENT DRIVE FOR FOUR QUBITS 


In this section we derive the driven cavity evolution for 
four qubits, as shown in FIG. 2 of the main text. The 
full system Hamiltonian is block diagonal, with one block 
for each state of the qubits. In each block the effective 
cavity Hamiltonian is 

He =Cbca^a + ao[cos (wdi^ + </'i) 

-I-cos (wd 2 ^ + </'2)]0(to - i) (a-I-a'^) , (6) 


where the shifted cavity frequency wc = <^c + XQ de¬ 
pends on the qubit-state dependent dispersive shift XQ = 
Xq((Tz), which defines the state of the qubits, and there¬ 
fore, the corresponding block of the full Hamiltonian. 
XQ = ±2 xq for odd-parity states and XQ = 0,±4 xq 
for even parity. 

We go to a frame rotating with the shifted cavity fre¬ 
quency (he to obtain 


= y©(^D - i) [a (e + e -k h.c.] , (7) 

where = £hc~i^Di- We formally solve for the evolution 
operator in each block for times t > to as 


U{t, 0) = Texp 





( 8 ) 


where T is the time ordering operator. The evolution 
operator for the full system then has the form 


3 


where Uj {t, 0 ) are the solutions to equation ( 8 ) for each 
qubit state \j) in the computational basis. 

We evaluate the integral in equation ( 8 ) using the Mag¬ 
nus expansion [34] to take care of the time ordering. Be¬ 
cause [a, a^] = 1 the Magnus expansion truncates at sec¬ 
ond order, and conveniently the second order term is a 
global phase that we can ignore. We evaluate the integral 
in two regimes, corresponding to even or odd parity, and 
in so doing obtain the solution for the blocks of the full 
Hamiltonian. For even parity, Ai ^ 0 and A 2 7 ^ 0, and 
to first order we have 







1 ) 


+ ^ - 1 ) ) + h.c. 


The evolution operator of equation ( 8 ) for even parity is 
then 


= b(a'E,{t-D)), ( 10 ) 


where 

aE(to) 


2 


giAiio _ 2 giA2tD _ \ 


Ai 


( 11 ) 


In general as (to) varies in phase for different even-parity 
qubit states as Ai _2 depend on the specific state. 

For an odd-parity qubit basis state, the shifted cavity 
will be on resonant with one of the two applied drives, 
and therefore, one A^ is zero, while the other is nonzero. 
Since the solution is symmetric in regards to which A^ = 
0 , we will examine Ai = 0 and A 2 7 ^ 0 without loss of 
generality. The integral in equation ( 8 ) now gives 



H'c{t')dt' 


2 


d 


^to + i 


g —■iA2tD 
^2 




which leads to the odd-parity evolution operator 

Uo{t, 0) = I)(ao(to)) = b (/?(to) + »7(to)), (12) 


where 

/3(to) = --yto (13) 

^(to) = -^(e*^^‘°-l). (14) 

In this case it is important to note that while ??(tD) varies 
between odd-parity states, /3(to) is the same for all odd- 
parity states. 

If wdi = 2xq and = —2xq as in the main text, 
then for even parity we have Ai,A 2 € {±2xq,±6xq}. 
For odd parity, we have chosen Ai =0, which means 
that A 2 = 4xq. Thus, if we set tu = tt/xq, we have that 
aE(to) = 0, and rj^tu) = 0, while /?(to) 0. As a result 

of this. 


(9) 


UEit,0) = b{0) = 1 
Uoit,0) = bmu)), 


(15) 

(16) 
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and the evolution operators are the same for all states 
within the same parity subspace. This, along with the 
block diagonal form of the full system Hamiltonian im¬ 
plies that the unitary on the full system will have the 
form 

Ud= 1-^'^ ® 0 ) + ^ \j) {j \0 Uo{t, 0 ) 

jGeven jGodd 

= P^^lc + P§^DiP), (17) 

where are projectors onto the even- and odd-parity 
subspaces respectively. 


We solve the master equation 

p{t) = - i[Hjc, p{t)] 

+ (^L^prnl - p(t)}) , (21) 

with 7 j = 200 MHz, 7 ^ = 200 MHz, 7 d = 1 MHz, 
and 5j/27r = 50 MHz, for both ojc ~ = 2xQ) and 

—2xq (where XQ/27r = 5 MHz) to generate the detection 
probability curves of the main text FIG. 3. 


MASTER EQUATION FOR THE CAVITY-JPM 
COUPLED SYSTEM 

In this section we describe the methodology of the nu¬ 
merical simulations used to obtain FIG. 3 of the main 
text, which shows the detection probability and parity 
contrast for our example threshold detector, the JPM. 

We start with the cavity-JPM coupled system for four 
qubit odd parity, where the shifted cavity frequency Cjq = 
wc ± 2 xq and the JPM frequency wj = wq are such that 
the JPM is ±2 xq detuned from the cavity. This leads to 
the Hamiltonian 

Hjc = wca'^'a - yUj -k gj (adl -k ) , 

where af couple the ground and excited state of the 
JPM. As shown previously [21, 28, 29, 35], the JPM 
is well-approximated as a three-level system with self- 
Hamiltonian Hj = —^(jj, where ctj = diag(l, —1, fc). 
The energy of the third ‘measurement’ state is arbitrary 
in our model; the JPM only tunnels into it incoherently 
at fixed rates [28, 29[. 

In addition, the JPM experiences a number of incoher¬ 
ent processes. Tunneling from the JPM excited state to 
the measured state corresponds to photon detection and 
occurs at the bright count rate 7 j, with the corresponding 
Lindblad operator 

^2 = v^(IIcG |m) (Ijj). (18) 

Inelastic relaxation takes the JPM from the excited state 
to the ground state at a rate 7 r, with corresponding 
Lindblad operator 

= v^(Ic G |0) (Ijj). (19) 

Finally, false detections, where the JPM tunnels from the 
ground state to the measured state, can occur at the dark 
count rate 7 d, with Lindblad operator 

^0 = v^(Ic G |w) (Ojj). (20) 

We ignore the effects of pure dephasing on the JPM as 
they do not affect the parity measurement protocol of the 
main text. 


QUBIT DISPERSIVE SHIFT MISMATCH 

In this section we examine in detail the possible sources 
of error caused by XQ mismatch that were described in 
the main text and shown in FIG. 4. Firstly, we quantify 
the decrease in measurement contrast and the coherence 
loss for two qubits. As described in the main text, this 
coherence loss can be restored by perfect reset. We then 
examine these effects for four qubits, and also quantify 
the change in the measurement basis due to detection 
probability mismatch, which is unique to four qubits. 


Two Qubits 


If the dispersive shifts for the two qubits vary, such 
that XQi = XQ and xQs = XQ + e, then the cavity-drive 
detunings within a parity subspace split, such that we 
now have 


Aqo = - 2 xq - e 
All = 2 xq + e 
Aoi = e 

Aio = -e ( 22 ) 

where A^- is the cavity-drive detuning for the state \ij). 

After a time = "^/XQ. the cavity will be in a qubit- 
state dependent coherent state with amplitude 


aoo = 

ail = 


ao 


2 ( 2 xq + e) 

gp / 
2 ( 2 xq + e) ^ 


(e ^ 


I —^ e 

e XQ - 1 


Q-o / ^ , 

aoi = (^e - 1 ) , 

gp / 1 

aio = -^ XQ _ 1 


(23) 


For the even-parity qubit states we have that an = 
—doo, and for the odd-parity states aio = — doi. As 
a result, states within the same parity subspace will 
have the same detection probability, and so there will be 
no intra-subspace decoherence during the measurement 
stage. However, as jaooj ^ 0 there will be a reduction 



of measurement contrast due to the increased probabil¬ 
ity of misidentification. As the coherent states for each 
parity subspace are out of phase with one another there 
also will be intra-subspace decoherence during the drive 
stage. We will now quantify these effects in the e/xQ < 1 
regime. 

For the measurement contrast and misidentification we 
can quantify the effect of xq mismatch by looking at the 
cavity occupation for odd parity, |aoP = |c«oiP = |cnoPj 
and for even parity, |q;eP = |aooP = |aiiP- For odd 
parity 



where e = e/xQ) we see that the reduction of 

|aoP from that when the dispersive shifts are perfectly 
matched is second order in e/xQ- Similarly, for even par¬ 
ity 



and the increase in |q;eP from that when the disper¬ 
sive shifts are perfectly matched is again second order 
in e/xQ- The fact that the lowest order dependence is 
quadratic originates from the expansion of the cosine, 
and can be physically explained by the fact that the am¬ 
plitude is originally tuned to zero, and a square-law de¬ 
tector responds to intensity which is |ap. 

To quantify the intra-subspace decoherence during the 
drive stage, we consider an arbitrary superposition of 
odd-parity states with the cavity initially in vacuum, to 
which we apply the modified qubit state dependent drive 
(accounting for XQ mismatch). The resulting state 

|T) =a|01)|aoi) + &|10)|aio), (26) 

where |ap -I- |6p = 1, is no longer a product state and as 
such there will necessarily be decoherence of the reduced 


qubit state. The reduced 

qubit state is 



0 

1 

0 

1 

0 


0 

ap Dab 0 

Pq = 

0 

Dab |6|2 0 


^0 

0 0 0 / 


with the complex decoherence factor 


D = (aoil aio) = exp {- (|aoiP + aoi)}- (28) 


From equation (23) we see that 


|aoi|^ = 2Ao ( 1 - cos ( —^ 

xqJ 


al, = A} 


7r6 \ 

1 — 2 ( cos ( — — i sin 


7re \ 


-I- cos 2— — isin 2 — 


XQj 


XQj 


XQJ \ 



( 29 ) 


(30) 


where Aq = aQf2€. 

Keeping only the first nonzero term in the Taylor series 
expansion we find 



As can be seen, for an odd-parity superposition, xq mis¬ 
match causes decoherence at a rate that is second order in 
the small parameter e/x, and also introduces a complex 
phase factor that is first order in e/x. 

For an even-parity superposition, we simply replace Aq 
with the corresponding expression for even states, and we 
find 


laool' 


Oq 


^00 


(2(2x0 





Thus the even parity situation is even better, as the deco¬ 
herence rate is now fourth order in e/x, and the complex 
phase factor now third order. 

For both odd- and even-parity superpositions full intra¬ 
subspace coherence will be returned to the two qubits by 
perfect reset, as it disentangles the qubits and the cavity 
in a unitary and phase insensitive way. However, if cavity 
decay occurs before reset, or imperfect reset leaves resid¬ 
ual photons which subsequently decay, the qubit state 
will lose coherence, as will be discussed in section . 


Four Qubits 

The situation is considerably more complex in the 
four qubit case as the cavity-drive detunings split into 
sixteen distinct frequencies, one for each qubit state. 
Correspondingly, after the qubit state dependent drive 
there are sixteen possible cavity states \aijki), where 
G {0,1} index the qubit state jijkl}. As a re¬ 
sult of the double frequency drive each will have 

two components oscillating at different frequencies, with 
each component similar in form to those of equation (23). 
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Cavity states for even qubit parity will have both com¬ 
ponents of the form of ago, where the amplitude is sup¬ 
pressed by a factor 1 /( 2 xq -|- e), as both drives are off 
resonance by at least 2xq ■ For odd qubit parity the cav¬ 
ity state will have one component that is similar in form 
to ago from the off resonance drive and one component 
similar to ogi from the nearly on resonance drive. 

Given the similar structure of aiju in the four qubit 
case to aij of the two qubit case it is clear that the errors 
introduced by xq mismatch in the four qubit case will be 
a generalization of those in the two qubit case. There will 
be a reduction of measurement contrast due to increase 
in the probability of misidentification, which will again 
be at most second order in the small parameters Ci/xQ- 
The overlap between any two states \aijki) corresponding 
to qubits in the same parity subspace is very similar to 
equations (31) and (32) for the two qubit case, and as 
a result, the intra-subspace decoherence during the drive 
phase will be at most second order in ti/xQ for odd- 
parity superpositions and fourth order for even parity. 
As in the two qubit case full coherence will be returned 
to the qubits by perfect cavity reset. 

Despite the XQ mismatch the two qubit case exhibits 
symmetry in within a parity subspace that ensures 
the detection probabilities are the same for cavity states 
corresponding to qubit states of the same parity. How¬ 
ever, in the four qubit case this symmetry is no longer 
present, as in general £2 7 ^ eg yf £ 4 , and as a result, de¬ 
tection probabilities will differ within a parity subspace. 

For example, consider an equal superposition of odd- 
parity states I a) and \b), which after the drive stage are 
entangled to cavity stages jaa/b) respectively. Detection 
by the photon counter results in back action on the cavity 
described by the non-unitary back action operator B. 
Defining the normalized states 


IV'a.b) = 




v£^a,b| B^ B |o;a,b) 


B I O-a^h) 

Pnh '' 


(33) 


where Pa,b is the detection probability of the state \aa,b), 
then the resulting state of the cavity-qubits system after 
detection by the photon counter is 

m = ^[VK\a)\i>a) + VPb\b)Hb)), (34) 

where the normalization factor N = Pa + Pb- 

The four qubit states of equation (34) are entangled to 
distinguishable cavity states such that the reduced four 
qubit density matrix’s coherence is reduced by the over¬ 
lap {il^b\ V’a)- To remove this effect and focus solely on 
detection probability mismatch we assume that we can 
perform perfect reset of the cavity and create the state 

!'!') = ^(/^l«) + ynl&))|0), (35) 


state |4 >)q is no longer the input state and this change is 
a result of the fact that different detection probabilities 
for states within the same parity subspace change the 
basis of measurement of the protocol. To quantify this 
effect we calculate the magnitude of the overlap between 
the target state |4't)q = (|a) -f \b)) and the state 
|4')q, given by 

0 = = + (36) 

which is unity for Pa = Pt as expected. 

If we assume the photon counter is a JPM with a sub¬ 
traction operator back action [30], then the detection 
probability of the state \ai) is 

Pi = 1 - exp {-[ofip} . (37) 


By Taylor expanding the overlap and discarding terms 
smaller than exp (— 2 |a£,/bP) we obtain the approximate 
overlap 

_ +!“'’! ) (cosh (|aaP - |abP) - l) 

-t O . (38) 

Since we are interested in the effect of detection prob¬ 
ability mismatch we can assume Ua/b differ in magnitude 
only and set \aa\ = |q;| and \ab\ = |a|(l-l-(I) without loss 
of generality. The overlap can then be rewritten as 

„ , e-l“l’("+“+‘“)(cosh(|a| 72 i + i2))-l) 

4(2-e-l“l"-e-l“l"(i+'5)") 


To examine the scaling of the overlap, we assume that 
d ^ 1 and only keep terms up to order <5^ in the Taylor 
series, which gives 


0 = 1 - 


4(l-e-l“l=) 


+ O (S^) . 


(40) 


Four qubit XQ mismatch causes a difference in the cavity 
coherent state amplitude for states in the same parity 
subspace that scales as oc (ef — e|)/xQ- As we assume 
ei/xQ is a small parameter, then (ef — e^)/XQ i® as small 
or smaller, so the assumption that 5 ^ 1 is valid. 

From equation (40) we see that the error in the output 
state scales as (ef — e|)/XQ and is damped by the fac¬ 
tor |a|^e“^l“l^/4. For |a| = 3 as used in the main text 
this damping factor is on the order of 10 “®, and so any 
error in the output state caused by detection probability 
mismatch is negligible. 

A better approximation to the detection probability in 
the case where the JPM can also relax inelastically is [5] 


for which the reduced qubit state is the pure state 
|'1')q = (VT^\a) + y/Pb\b)) /VN■ Clearly, the output 




1 — exp 


7J 1 
7J + 7R J ' 


(41) 
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Using 7 j = 7 r as elsewhere, then the error of equation 
(40) is damped by a factor oc . For |q;| = 3 as 

before, this is on the order of 10 “^, which again results in 
a negligible effect that will only decrease with increasing 
cavity photon number. 

QUBIT DECOHERENCE AND CAVITY DECAY 

When xq mismatch is present there can be a phase 
difference between cavity states corresponding to qubit 
states with the same parity, as discussed for two qubits in 
section . In this section we examine the effect this has on 
intra-subspace coherence when a cavity decay mechanism 
is introduced, as described in the main text. 

We start with a typical post drive or imperfect reset 
stage state of the form 

P(0) = ® Ps.s'(O) ® Pb, (42) 

s,s' 

where s,s' label the state of the qubits, /5s_s'(0) are the 
qubit dependent cavity matrices, and the environmental 
bath is represented by the state pb, which is assumed to 
be stationary and a thermal state. Additionally we as¬ 
sume that the qubit states have no intrinsic time depen¬ 
dence, i.e. we work in the rotating frame of the qubits. 

We consider the intial state of the qubits to be an equal 
superposition of two states in the same parity subspace, 
which for simplicity we have labelled as |0) and |1). This 
results in the initial system state 

Ps(0) = ^( |0) (0| ® l“o) (aol + |0) (1| 0 |ao) (ai| 

-f |1) (0| (g) |ai) (aol + |1) (1| <8 |ai) (ai| ), (43) 

where Ps,s'( 0 ) = (c^s'l- Due to xq mismatch ag and 

ai differ in their phase only and we can set 

ao = |Q:|e*‘^“, ai = (44) 

Assuming that the cavity is in a coherent state is accurate 
for a post drive stage state, but the post reset state will 
be a more general state of the form ps^s'{Q) = IV's) (V's'l • 
Nevertheless, we are eventually interested in steady state 
dynamics which are independent of the initial state of 
the cavity (under the Born-Markov approximation). We 
will treat the post drive stage situation first and then 
generalize to the post reset state. 

To solve for ps^s’it), we will use the Wigner character¬ 
istic function approach. We define the Wigner character¬ 
istic function Xss'{c(,t) of the state Ps,s'{t) as 

Ps,s'{t) = ^ J (faXss'{a,t)D{-a), (45) 

where D{P) is the displacement operator. 


If we assume a Jaynes-Cumming type interaction be¬ 
tween the cavity and a bosonic environment with smooth 
spectral density, then following [36] we arrive at the an¬ 
alytic solution for Xss' (o^j t) 

Xs«'(a,t) =X..'(«e-‘("-*“),0)e(3l“l=U-=*'‘-i)), 

where k is the decay rate of the cavity, w is the frequency 
of the cavity mode, and p = l + 2n{uj) with n(w) the Bose 
distribution at frequency w. 

In order to quantify the absolute maximum amount of 
decoherence, we are interested in the steady state solu¬ 
tion (t —>• oo) of equation (46) 


Xas'(a) = Xss'(0,0)e 

(47) 

Using the facts that 


Xsa'(0,0) = y;Tr[ps_^/(0)] 

(48) 

Xio{a,t) = xoi{-a,t) 

(49) 

we obtain 


Xoo(a) = Xii(a) = 

Att 

(50) 

Xoi(a) = exp - 1)) 

(51) 

Xio(a) = exp - 1)^ 

(52) 

Each of equations (50), (51) and (52) can 

duced to Xss' {a) = (oq, ai)XThermal(a, w). 

be re- 
where 


XThermai(a,a;) = e“2l“l^/47r is the Wigner characteris¬ 
tic function for a thermal state at frequency w. 

In light of this simplification the steady state solution 
for the full system reduces to 

PS = '^ c(s, s')Fss' {ao, ai) |s) (s| 

s,s' 

® “ y (i^Q;XThermai(a,a;)ZI(-a), (53) 

where = exp — 1)). As this 

is a product state it is the trivial to trace out the state 
of the cavity and obtain the state of the qubits alone 

PQ = l{\0) (0|+Foi(ao,ai)|0)(l| 

+ Fio(ai,ao)|l)(0| + |l)(l|). (54) 

For two qubits, we examine the damping envelope of 
the off diagonal elements, given by 

Foi{ao,ai) = exp (-|ap(e^*‘^“ -f 1)), (55) 

where we have used the fact that (po — pi = 2ipo ~ 
for two qubits, as can be calculated from either pair of 
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expressions in equation (23). We are interested in the 
absolute value of this envelope, which is given by 


|Foi(Q!o,ai)| = exp (-lapRe + 1)}) 


= exp 




(56) 


where the second equality comes from using equation (23) 
to define . From this we see that the decoherence de¬ 
pends on both the cavity occupation and the magnitude 
of the xq mismatch, being maximal when e/xq = 1. Us¬ 
ing |ap = 9 as used elsewhere, and a strong mismatch of 
e/xq = 0.1, the qubit coherence is reduced to 64% in the 
steady state. However, for a high-Q cavity we do not ex¬ 
pect any photon loss during the parity readout protocol, 
and as such, the post reset state is of greater interest. 

Generalizing to the post reset state we need only 
change the definition of the damping envelope to 


Fss' =TT[ps,s'm = {ilJs'\^s), (57) 

where \ips) is the state after detection and imperfect reset 
for an initial cavity state |as). We will consider the worst 
case scenario for the states |q;o/i)j where the damping 
envelope of equation (56) reduces to exp(—2|ap). This 
occurs for e = XQj where the states |ao/i) li^ on the real 
axis with ag = —ai. This is the worst case as the states 
are maximally distinguishable, and correspondingly we 
expect the states |V'o/i) to also be maximally distinguish¬ 
able in this case. 

Figure 5 shows a numerical simulation of 1 — Foi for 
the worst case post reset state as a function of both the 
initial cavity occupation jap, and the number of photons 
removed by the photon detector. For these simulations, 
the detector was assumed to be a JPM with subtraction 
operator back action [30]. Multiple photons can be re¬ 
moved due to JPM inelastic relaxation. As can be seen, 
even for the worst case scenario, qubit coherence is re¬ 
duced by less than 1%. As we expect e ^ XQ ^ realistic 
experiment, the decoherence caused by post reset state 
photon loss will be far from the upper bound presented 
here, and can be considered inconsequential. 


HIGHER ORDER EFFECTS BEYOND THE 
DISPERSIVE HAMILTONIAN 

While the dispersive Hamiltonian is the first order ap¬ 
proximation to the cavity-qubit coupling in the relevant 
regime of our protocol, it is worthwhile to consider the ef¬ 
fects on the protocol of the full Jaynes-Cummings Hamil¬ 
tonian for the cavity-qubit coupling, given by 

N ^ N 

Hjc = wco'^'a-^ + (d+a -f , (58) 

fc =0 fc =0 



FIG. 5: 1 — fbi as a function of the cavity occupation prior 
to detection and imperfect reset (horizontal axis), and the 
number of photons removed from the cavity by the photon 
detector (color and line style). The initial cavity states con¬ 
sidered are maximally distinguishable, and so this represents 
the worst case scenario. 


where is the cavity-qubit coupling for qubit k. The 
major effects of the full Jaynes-Cummings Hamiltonian 
occur during the drive stage. 

The first is an asymmetric shift to the cavity frequency 
within a parity subspace, due to the Kerr-like interaction 
term, which is one order of approximation higher than 
the dispersive Hamiltonian. The result of this is that the 
minima of the even qubit parity cavity occupation curves 
no longer coincide, as can be seen in FIG. 6(a) for the 
four qubit case. The second major effect is the formation 
of dressed cavity-qubit states, the result of which is that 
there is residual cavity occupation for even-parity states, 
and the minima are no longer zero for all even-parity 
states. Therefore, for no choice of to will the cavity be in 
the vacuum state if the qubits are in an even-parity state. 
This reduces the measurement contrast, though as can be 
seen in figure FIG. 6(b) this reduction is only on the order 
of a few percent, and measurement contrasts approaching 
92% can still be achieved with the same parameters as in 
the main text. 

These effects can be mitigated by increasing the cavity- 
qubit detuning, in which case the parity measurement 
contrast asymptotes to the dispersive value. However, 
this either decreases XQ ^-nd therefore increases the 
length of the drive stage, or necessitates also increasing 
the cavity-qubit couplings (to keep xq constant), which 
is only advisable to the point where the rotating wave ap¬ 
proximation begins to break down. It is also possible to 
improve the contrast by better pulse shaping in the drive 
stage, and this will be the focus of future work. As can 
also be seen in FIG. 6(a), the fact that the even qubit 
parity cavity photon states are not the same between 
parity bands would cause intra-subspace qubit decoher¬ 
ence. Full coherence would be returned by perfect reset, 
unless photon loss occurs (either through a detection or 
other mechanism). However, photon loss of any kind is 
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so unlikely for these low occupation cavity states that the 
effect on qubit coherence is minimal, on the order of a 
few percent at most. 


FIG. 6: (a) Cavity occupation for even-parity qubit states 
with Jaynes-Cummings cavity-qubit coupling, (b) Bright and 
dark count rates, and measurement contrast for the cavity 
photon numbers of (a), assuming the even state is from the 
right band (worst case). 













